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I. 


INTRODUCTION 


In  recent  years,  tremendous  advances  have  been  made  in  the  basic 
understanding  of  shaped  charge  mechanics.  These  include  the  processes 
of  liner  collapse,  jet  formation,  and  formulas  for  jet  radius  and 
strain.  Both  two-dimensional  computer  codes  and  one-dimensional  model- 
ing have  been  used  successfully  for  the  analysis  of  shaped  charges. 

One  important  area,  however,  has  not  been  well  understood  and  cannot 
be  simulated  easily  by  computer  methods.  This  area  is  the  breakup  of 
the  shaped  charge  jet.  This  breakup,  or  segmentation,  occurs  in  all 
Jets  that  have  a large  velocity  gradient  with  higher  speeds  at  the  tip 
and  lower  speeds  near  the  tail.  This  process  can  be  seen  in  Fig.  1 
which  shows  flash  radiographs  of  a typical  shaped  charge  jet  at  three 
successive  times,  displayed  in  the  proper  position-time  coordinates. 

At  the  earliest  time  the  jet  segmentation  has  begun  near  the  tip;  most 
of  the  Jet  is  still  continuous.  The  later  times  show  the  completely 
segmented  jet. 

In  [1]  and  [2],  we  have  presented  initial  studies  on  the  jet 
breakup  phenomenon.  Let  us  now  briefly  review  the  results  of  these 
Initial  studies.  Formulas  for  the  strain  and  radius  of  shaped  charge 
jets  based  on  a one-dimensional  model  are  presented  in  [1],  Further, 

in  [1]  and  [2],  a method  to  determine  the  jet  breakup  time  distribu- 
tion from  timed  flash  jet  radiographs  was  developed.  This  method  was 

then  applied  to  jets  from  a series  of  identical  BRL  81.3mm  standard 

42°  copper-lined  charges  and  an  unconfined  105mm  42°  copper  lined 
charge  in  [1].  Despite  the  scatter  in  the  results,  the  breakup  time 
does  show  a definite  "trend"  as  is  indicated  in  [1],  In  fact,  for  the 
particular  charge  studied,  the  trend  indicates  that  the  jet  breaks 
first  near  the  tip  with  a progressively  increasing  breakup  time  towards 
the  tail.  Thus,  shaped  charge  jets  do  not  necessarily  break  simulta- 
neously along  their  length,  as  assumed  previously  by  some  investigators. 
The  resulting  jet  breakup  time  distribution  was  then  contrasted  to  the 
one-dimensional  theoretical  jet  strain  and  radius  distributions.  This 
is  presented  in  [1]  and  [2]  where  it  is  concluded  that,  for  copper 
liners,  the  breakup  tiiie  distribution  is  related  to  the  jet  radius  dis- 
tribution. No  other  correlation  was  indicated  in  [1]  or  [2]. 

The  present  study  addresses  the  problem  of  jet  breakup  through  the 
use  of  two  approaches:  (i)  triple-flash  radiographs  of  various  jets  to 
determine  jet  velocity,  breakup  time,  etc.,  and  (2)  surface  instability 
as  a cause  of  breakup. 

The  first  approach  is  actually  an  extension  or  continuation  of  the 
work  we  have  presented  in  [1].  The  method  of  [1]  to  measure  breakup 
time,  jet  velocity,  etc.  from  timed  flash  radiographs  is  applied  here 
to  many  jets  from  various  shaped  charges  having  different  cone  angles, 
different  liner  wall  thicknesses,  liners  with  tapered  walls,  different 
liner  materials,  and  light  or  heavy  confinement,  This  approach  has 
enabled  us  to  obtain  many  characteristics  of  jets  from  a large  cross- 
section  of  shaped  charges.  The  data  from  all  of  these  charges  are 
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conveniently  tabulated  and  graphically  displayed  in  this  report. 

The  geometry  of  the  chargee  studied  are  given  in  Appendix  A.  The 
measured  jet  velocity,  breakup  time  and  jet  radius  are  plotted 
graphically  in  Appendix  B.  Appendix  C contains  a tabular  list  of 
this  data  and  in  addition  includes  the  aspect  ratio  of  the  jet  seg- 
ments and  difference  in  velocity  between  neighboring  particles.  All 
of  this  data  has  been  analyzed  to  obtain  trends  in  the  breakup  phe- 
nomenon. The  breakup  trends  have  been  contrasted  with  the  one- 
dimensional  strain  and  radius  predictions  using  the  formulas  of  [1]. 
Within  an  individual  jet  no  general  correlation  has  been  observed 
between  breakup  time  and  one-dimensional  parameters.  The  results 
of  the  breakup  time  measurements  of  all  copper  jets  studied,  however, 
do  exhibit  a definite  trend.  This  trend 

yields  a semi-empirical  design  curve  for  breakup  time  which  can  be 
used  in  conjunction  with  a one-dimensional  shaped  charge  model  to 
predict  breakup  time  for  copper  jet9.  This  first  approach  is  pre- 
sented in  Section  II. 

The  second  approach  involves  the  study  of  the  instability 
caused  by  various  disturbances  in  the  jet  with  the  goal  of  deter- 
mining whether  or  not  this  instability  may  cause  the  jet  to  breakup. 

We  know  from  classical  hydrodynamics  that  a continuous  liquid  jet 
will  break  into  small  segments  because  of  surface  instability.  The 
breakup  of  a shaped  charge  jet  resembles  this  very  much,  and  our 
study  shows  that  the  shaped  charge  jet  is  indeed  subject  to  surface 
instability.  In  principle,  there  are  three  forces  which  may  cause 
this  instability:  surface  tension,  aerodynamic  force,  and  material 
strength  (elastic-plastic  force)  in  the  jet.  Both  analytical 
studies  and  two-dimensional  finite-difference  numerical  calculations 
were  made,  and  the  results  indicate  that  material  strength  is  the 
main  cause  for  shaped  charge  jet  instability  and  breakup.  In  addition, 
the  effect  of  strain-rate,  time  of  disturbance  initiation,  and  inertia 
were  considered.  The  results  of  this  stability  study  are  discussed 
in  Section  III.  Finally,  general  conclusions  are  stated  in  Section  IV. 

We  would  like  to  note  that  the  results  presented  in  this  report 
summarize  the  work  conducted  during  the  entire  contract  period.  Some 
of  these  results  were  presented  previously  in  quarterly  progress  re- 
ports. 
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II.  JET  RADIOGRAPH  DATA 


The  results  which  we  presented  in  [1]  indicate  that  the  measure- 
ment of  Jet  radiographs  provides  an  efficient  and  useful  procedure  to 
obtain  breakup  data.  Also  breakup  trends  are  indicated  in  [1]  which 
warrant  further  investigation.  In  this  study,  we  have  theiefore  con- 
tinued using  this  approach  to  study  jets  from  a large  number  of 
charges  having  different  geometries  and  liner  materials. 

A.  General  Approach 

The  detailed  development  and  equations  for  the  computation  of 
breakup  time.  Jet  velocity,  etc.  from  the  radiograph  measurements  are 
given  in  [1].  Here,  we  will  give  only  a brief  description  of  this 
method.  Figure  2 shews  the  position-time  plot  of  two  typical  neigh- 
boring Jet  segments  labeled  as  K and  K+l.  Suppose  we  have  radiographs 
of  these  two  segments  at  times  tj^  and  C2  a®  indicated  in  Figure  2. 
(Note:  in  the  actual  case,  we  have  analyzed  three  timed  flash  radio- 
graphs, as  supplied  by  BRL) . We  may  then  compute  the  velocity  of  each 
segment  by  measuring  how  far  the  segments  have  traveled  during  the 
time  t2~tp  In  the  two  radiographs  the  segments  are  separated  by  a 
gap  as  shown.  Under  the  assumption  that  the  segments  remain  at  con- 
stant length  and  constant  velocity  after  breakup  (which  is  a good 
approximation  after  the  jet  is  completely  broken  as  evidenced  by  the 
radiographs),  it  is  a simple  matter  to  trace  back  the  front  of  seg- 
ment (K+l)  and  the  rear  of  segment  K until  they  meet,  i.e.  the  gap  be- 
comes zero.  The  time  of  this  meeting  t^  is  the  breakup  time  of  seg- 
ment K from  segment  (K+l).  This  procedure  may  be  applied  successively 
to  each  pair  of  particles,  i.e.  K and  (K+l),  (K+l)  and  (K+2),  etc., 
until  the  complete  breakup  time  distribution  is  obtained.  There  is  a 
certain  degree  of  approximation  to  this  procedure.  Since  the  particles 
do  not  break  successively  from  one  end  of  the  jet  to  the  other,  we 
sometimes  have  the  situation  where  a larger  segment  breaks  from  the 
main  jet,  continues  to  stretch,  and  then  finally  breaks  into  smaller, 
constant  velocity,  constant  length  segments.  However,  the  time  between 
the  initial  break  and  the  succeeding  breaks  is  small,  therefore  the 
procedure  may  still  be  applied  approximately  to  this  breakup  situation. 
The  resulting  breakup  times  have  been  found  to  be  within  the  accuracy 
of  the  experimental  procedures. 

This  method  of  determining  the  breakup  time  was  then  applied  to 
jet  radiographs  from  many  various  shaped  charges.  These  radiographs 
were  obtained  from  BRL  through  the  courtesy  of  R.  Jameson,  R.  Karpp, 

J.  Simon  and  J.  Majcrus.  For  convenience,  the  different  charges  are 
listed  In  Table  I;  the  detailed  geometry  and  drawings  ot  each  are 
given  in  Appendix  A.  It  can  be  seen  from  this  table  that  we  have 
studied  the  effects  of  cone  angle,  wall  thickness,  wall  taper,  and 
liner  material  on  the  breakup  mechanism. 


Table  1 


List  of  Chargee  Examined* 


Copper  Lined 

Charge 

Cone 

Diameter 

Wall  Thickness 

Mo. 

Angle 

mm 

sa 

1 

20* 

38.1 

1.168 

2 

40* 

38.1 

1.168 

3 

60“ 

38.1 

1.168 

4 

90“ 

38.1 

1.168 

5 

42* 

50.8 

0.762 

6 

42“ 

50.8 

1.524 

7 

42* 

50.8 

2.540 

8 

42*** 

81.3 

tapered  wall 

Aluminum  Lined 

Charge 

Cone 

Diameter 

Wall  Thickness 

No. 

Angle 

mm 

mm 

9 

40* 

38.1 

1.626 

10 

60* 

38.1 

1.626 

11 

90* 

38.1 

1.626 

* Note  that  we  have  also  performed  a breakup  study  of  jets  from  the  BRL 
Standard  81.3ma  charge  and  the  105mm  unconfined  charge  In  [1]. 

**  Angle  of  inside  wall 
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B.  Results  of  Radiograph  Measurements 

From  the  measurements  of  the  jet  radiographs,  we  have  determined 
the  jet  velocity,  jet  breakup  time,  and  jet  radius  distributions  for 
all  of  the  charges  listed  in  Table  I.  Plots  of  all  of  the  data  ob- 
tained are  given  in  Appendix  B.  Figure  31  shows  the  jet  velocity 
distribution  along  the  length  of  the  jet  at  a particular  time  for  the 
series  of  31.8mm  copper  lined  charges  having  cone  angles  of  20s,  40°, 
60s  and  90s  (Charge  Nos.  1-4).  Test  results  of  these  charges  were 
first  reported  in  [3],  where  the  average  breakup  time  of  the  complete 
jet  was  given.  The  time  base  used  in  all  of  the  present  studies  is 
given  as  t ■ 0 when  the  detonation  wave  reaches  the  apex  of  the  cone.* 
From  Fig.  B1  we  see  that,  for  a given  liner  material,  the  tip  velocity 
of  the  jet  decreases  as  the  cone  angle  increases.  It  should  be  noted 
that,  for  the  20s  copper  cone,  the  resulting  jet  possessed  a bifur- 
cated region  at  the  tip  portion.  We  therefore  made  our  radiograph 
measurements  starting  from  the  first  fully  coherent  particle  of  the 
jet,  which  has  a velocity  of  8.4mm/usec,  much  smaller  than  9.9mm/usec 
reported  in  [3].  We  should  further  note  that  it  has  been  determined 
that  this  bifurcated  region  is  caused  by  supersonic  flow  in  the  jet 
formation  region  (See  [4]  for  details).  None  of  the  other  charges 
examined  under  this  study  possessed  this  bifurcated  region. 

The  jet  velocity  distributions  for  the  remaining  charges  are 
shown  in  Figures  B2-B6.  It  is  interesting  to  note  that  Fig.  B6  in- 
dicates that  the  81.3mm  tapered  wall  charge  studied  here  has  a higher 
jet  velocity  (8.3  km/sec  at  the  tip)  than  the  81.3mm  standard  charge 
with  constant  wall  thickness  studied  in  [1]  (7.9  km/sec  at  the  tip). 
This  somewhat  surprising  result  has  a fairly  simple  explanation.  Be- 
cause of  the  wall  taper  (thicker  at  the  cone  apex  tapering  linearly 
toward  the  cone  base),  as  we  proceed  from  the  apex  toward  the  base  the 
liner  elements  reach  higher  collapse  velocity  than  corresponding  ele- 
ments of  the  charge  with  a constant  wall  thickness.  This  effect  pro- 
vides the  tapered  wall  liner  with  an  overall  decrease  in  the  collapse 
angle  (i.e.  the  angle  of  the  collapsing  liner  with  the  charge  axis, 
denoted  as  3 in  (1)).  From  the  basic  one-dimensional  jet  formation 
theory  [5]  a decrease  in  collapse  angle,  in  general,  produces  a 
higher  jet  velocity.  Finally,  from  iigs.  B1-B6  we  also  observe  that 
the  velocity  distribution  for  each  of  the  charges  studied  is  approxi- 
mately a straight  line. 

Plots  of  the  breakup  time  of  the  jet  particles  versus  their 
position  in  the  jet  at  a particular  time  for  the  eleven  charges  ex- 
amined are  given  in  Figs.  B7-B17.  The  scatter  in  the  results  Is  to 
be  expected  because  of  the  nature  of  the  breakup  mechanism.  In  all 
cases  a definite  trend  in  breakup  time  exists.  To  demonstrate  this 
trend,  a straight  line  was  passed  through  the  data  using  the  method 
of  least  squares.  This  line  is  also  shown  in  Figs.  B7-B17.  All  of 
the  jets  studied  exhibit  trends  which  indicate  that  the  tip  portion 
breaks  earlier  than  the  rear  portion  except  for  the  38.1mm,  40° 
copper  charge  (Charge  No. 2)  and  the  81.3mm,  42°  copper  tapered  wall 


* 


Except  for  the  data  which  is  tabulated  in  Appendix  C:  The  data  there 
is  based  on  the  initiation  time  of  the  main  explosive  charge.  See 
Appendix  C for  details. 


charge  (Charge  No. 8).  Breakup  trenda  for  these  two  charges  show  that 
these  jets  breakup  almost  simultaneously.  Note  that  since  only  part 
of  the  jet  for  each  charge  appeared  In  the  radiographs,  these  con- 
clusions are  only  valid  for  the  portion  of  the  jets  measured. 

Next,  the  measured  jet  radius  distributions  for  the  eleven  charges 
studied  are  given  in  Figures  B18-B28.  We  see  from  Figs.  B18-B21  that 
for  copper  lined  charges  of  a given  diameter  and  wall  thickness  the 
average  jet  radius  increases  with  increasing  cone  angle.  This  effect 
is  not  so  pronounced  in  the  case  of  the  aluminum  lined  charges  as  can 
be  seen  from  Figs.  B22-B24.  Further,  we  note  from  Figs.  B25-B27  that, 
all  else  being  equal,  relatively  large  increases  in  wall  thickness 
give  only  small  increases  in  jet  radius. 

Finally,  all  of  the  data  resulting  from  the  radiograph  measure- 
ments are  tabulated  in  Appendix  C.  In  addition  to  the  jet  velocity, 
breakup  time,  and  Jet  radius,  which  were  graphically  displayed  in 
Appendix  B,  we  have  also  tabulated  the  aspect  ratio  (length  to  diameter, 
i/d)  of  each  segment  and  the  jet  velocity  difference  between  neighbor- 
ing segments  AV,.  The  average  i/d  and  average  AVj  are  summarized  in 
Tables  II  and  III.  In  these  tables,  we  have  also  summarized  this 
data  for  the  iJRL  Standard  81.3mm  and  105mm  unconfined  charges  which 
were  studied  in  [1].  Since  there  was  a certain  amount  of  scatter  in 
these  quantities,  the  coefficient  of  variation  (standard  deviation  di- 
vided by  tha  mean)  in  each  case  was  also  computed.  From  Table  II  we 
see  that  the  average  i/d  varies  between  4 and  6 for  copper  lined 
conical  charges  of  constant  wall  thickness.  The  tapered  wall  cone 
has  an  i/d  ratio  of  3.3  which  may  indicate  that,  even  though  this 
charge  has  a higher  tip  velocity  than  the  corresponding  constant  wall 
charge,  its  penetration  performance  at  longer  standoffs  may  be  poorer 
because  of  the  small  segment  size  and  early  breakup  time.  The  alumi- 
num lined  charges  studied  have  average  i/d  ratios  of  3.89,  4.92,  and 
8.39. 


From  Table  III  we  observe  that  the  overall  average  velocity 
difference  between  neighboring  segments  for  the  copper  lined  charges 
Is  approximately  llOm/sec.  Recently,  Held  [6]  has  published  that  the 
average  AV,  for  the  copper  lined  German  charges  he  studied  was  approx- 
imately lOOm/sec.  For  the  aluminum  lined  charges  studied  here,  we 
found  an  overall  average  AVj  of  145m/sec.  Finally,  we  obtained  radio- 
graphs of  a copper  lined  French  charge,  the  ISL  "S2T",  from  Perez  [7]  and 
found  the  average  i/d  to  be  5.62  and  the  average  AVj  to  be  137m/sec.  Data 
of  this  nature  is  useful  in  obtaining  rough  estimates  of  the  number  of 
jet  segments  in  newly  designed  charges. 

C.  Comparison  of  Data  with  Theory 

In  [1]  and  [8]  an  improved  one-dimensional  theoretical  shaped 
charge  model  was  developed.  This  improved  model  is  based  on  the 
classical  Pugh-Eichelberger-Rostoker  [5]  theory  but  uses  a semi- 
empirical  collapse  formula  to  describe  the  explosive-metal  interaction 
process.  Further,  formulas  for  jet  strain  and  jet  radius  were 
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Table  II 


Copper  Lined 

Charge 

No. 

Average  Aspect  Ratio  (i/d)  of 

Jet  Segments 

Coeff.  f 
Variacion 

Geometry 

(Dia.,  Cone  Angle, 
Wall  Thk.) 

1/d 

(Average) 

1 

38.1  on,  20*,  1.168  am 

4.69 

.520 

2 

38.1  am,  40*,  1.168  am 

5.95 

.411 

3 

38.1  am,  60*,  1.168  ma 

5.91 

.447 

4 

38.1  ton,  90*,  1.168  am 

4.49 

.301 

5 

50.8  ma,  42*,  0.762  aa 

4.18 

.412 

6 

50.8  nun,  42*.  1.524  or. 

4.48 

.361 

7 

50.8  ma,  42*,  2.540  ma 

4.48 

.336 

8 

81.3  mm,  42*,  tapered 

3.31 

.417 

* 

81.3  ma,  42*.  1.905  mm 

5.93 

.375 

* 

86.4  mm,  42®,  2.921  mm 

5.21 

.324 

Average  i/d  for  all  copper  jets 

4.86 

c.v. 

.180 

Aluminum  Lined 


Charge 

No. 

Geometry 

(Dia.,  Cone  Angle, 
Wall  Thk.) 

i/d 

Coeff.  of 
Variation 

9 

38.1  mm,  40°,  1.626  mm 

3.89 

.454 

10 

38.1  mm,  60®,  1.626  mm 

4.92 

.420 

11 

38.1  mm,  90°,  1.626  mm 

8.39 

.426 

Average  i/d  for  all  Ai  jets 

5.  73 

C.V. 

.411 

* These  charges  were  studied  in  [1]  and  this  data  is  included  here 
for  completeness 
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Tabic  111 


Avaraga  Jat  Velocity  Difference  Between  Neighboring  Segments  (AVt) 


Copper  Lined 

Geometry 

Charge 

(Dla. . 

Cone  Angle, 

AV} 

Coeff.  of 

No. 

Wall  Thk. ) 

(m/sec) 

Variation 

1 

38.1  mm. 

20*, 

1.168  ma 

106 

.474 

2 

38.1  mm. 

40* , 

1.168  mm 

111 

.367 

3 

38.1  tma. 

60\ 

1.168  ma 

135 

.335 

4 

38. 1 am. 

90* , 

1. 168  ma 

109 

.344 

5 

50.8  mm. 

42*. 

0. 762  am 

110 

.497 

6 

50.8  mm. 

42*. 

1.524  mm 

115 

.522 

7 

50.8  mm. 

42*. 

2.540  ma 

119 

.449 

8 

81.3  mm. 

42\ 

tapered 

108 

.519 

* 

81.3  ma. 

42*. 

1.905  mm 

96 

.557 

* 

86.4  mm. 

42*. 

2.921  mm 

113 

.445 

Average 

AV  for  all  copper  jets 

112 

J 

C.V. 

.090 

Aluminum  Lined 

Geometry 

Charge 

(Dla., 

Cone  Angle, 

AVj 

Coeff.  of 

No. 

Wall  Thk.) 

(m/sec) 

Variation 

9 

38.1  ma. 

40*, 

1.626  sxa 

143 

.362 

10 

38.1  mm. 

60°, 

1.626  mm 

130 

.528 

11 

38.1  on. 

90°, 

1.626  mm 

162 

.154 

Average 

AV  for  all  A l jets 

145 

J 

C.V. 

.111 

* These  charges  were  studied  in  [1]  and  this  data  is  included  here  for 
completeness. 
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developed  for  this  model.  Since  ehe  publication  of  [1]  end  [8],  we 
have  Improved  this  model  further  by  Incorporating  the  effect  of  liner 
acceleration  during  collapse.  This  enables  the  prediction  of  the  In- 
verse jet  velocity  gradient  region  during  formation. 

This  improved  one-dimensional  model  was  applied  to  the  eleven 
charges  studied  here.  In  Figs.  B1-B6  the  theoretical  jet  velocity 
distribution  of  each  charge  is  shown  in  contrast  to  the  experimental 
data.  We  see  that  excellent  agreement  is  obtained  between  theory  and 
experiment.  Note  that  in  Fig.  B4  and  B5  there  is  some  discrepancy 
between  the  experimental  tip  velocity  and  the  tip  velocity  predicted 
by  the  theory  for  Charge  Nos.  6 and  7.  This  prediction  could  be  im- 
proved by  using  a different  value  for  the  acceleration.  At  present, 
the  acceleration  values  used  in  all  the  cases  are  given  by  a simple 
empirical  formula  which  depends  only  on  liner  density  and  thickness. 

For  these  two  particular  charges  the  liner  thickness  to  liner  diameter 
ratio  is  much  higher  than  the  other  chargee,  studied  and  therefore  thx.: 
simple  acceleration  formula  may  be  inaccurate  for  these  two  charges. 
Theoretical  jet  radius  distributions  were  also  computed  from  the  one- 
dimensional model  at  the  breakup  times  t^.  The  tj,  lines  in  Figs.  B7- 
B17  were  used.  As  shown  in  Figs.  B18-B28,  the  theoretical  radius  dis- 
tributions compare  reasonably  well  with  experimental  values. 

We  next  compared  the  breakup  time  distribution  with  the  theoret- 
ical distribution  of  other  jet  properties.  The  goal  here  was  to  find 
some  correlation  among  these  distributions  which  would  indicate  con- 
trolling parameters  in  the  breakup  mechanism.  We  compared  trends  in 
the  following  properties  to  the  breakup  trend:  let  strain  n,  jet  radius 
r j , and  jet  strain  rate  n.  The  quantities  n and  rj  are  defined  in 
[1],  and  n is  simply  the  first  time  derivative  of  n defined  there.  W'e 
also  contrasted  these  properties  to  the  amount  of  time  elapsed  from 
the  formation  of  an  element  until  it  breaks  up.  This  quantity  is  the 
absolute  breakup  time  minus  the  absolute  time  when  the  element  is  first 
formed,  and  is  denoted  t^-tf.  The  quantity  t^  is  from  the  experimental 
least  squares  line  and  the  quantity  t£  is  computed  from  the  one- 
dimensional model. 

Since  the  quantities  jet  strain,  jet  radius,  and  jet  strain  rate 
before  breakup  are  continuously  changing  with  time,  it  is  appropriate 
to  compare  the  breakup  time  of  each  segment  with  these  quantities  at 
the  time  when  each  segment  breaks  up.  To  do  this,  we  first  trace  back 
each  segment  in  the  radiograph  to  its  original  position  in  the  cone, 
x,  using  the  one-dimensional  model.  We  are  then  able  to  compute  n,  rj , 
and  ' of  the  segment  at  its  own  particular  breakup  time  from  the  least 
squares  experimental  trend.  Figure  3 shows  a plot  of  n,  rj , n,  tb  and 
tb~tf  versus  x/h  for  the  38.1mm,  60°  copper  charge  (Charge  No. 3).  Note 
that  x denotes  the  original  liner  position  and  h denotes  the  original 
height  of  the  cone.  Also  plotted  in  Fig.  3 is  rj  at  t * 55  isec  which 
is  a time  before  any  breakup  has  occurred. 


1 1 
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RELATIVE  CONE  POSITION  (x/h) 

FIGURE  3.  Theoretical  jet  strain-rate,  strain, 
radius  and  experimental  breakup  time 
vs.  relative  cone  position  for  the 
38.1mm,  ftO°  Cu  charge. (Charge  No.  3) 
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Figure  3 indicates  that  there  is  no  reasonable  correlation  be- 
tween the  breakup  time  distribution  and  n at  t^  or  n at  tb*  Both  the 
rj  at  t],  end  rj  at  S3  usee  appear  to  follow  the  general  trend  of  the 
breakup  time,  riots  similar  to  Figure  3 for  the  remaining  10  charges 
were  also  plotted  but  are  not  presented  here  in  the  interest  of  space*. 
However,  the  results  indicate  that  within  one  jet  even  r^  at  tb  and  rj 
at  a time  before  breakup  may  not  follow  the  trend  of  breakup  time. 
Results  do  indicate  that  for  copper  lined  charges  rj  at  a time  before 
breakup  seems  to  Increase  as  the  breakup  time  increases  but  no  general 
conclusion  may  be  drawn. 

We  next  looked  at  the  breakup  data  from  all  the  copper  lined 
charges  together  rather  than  just  throughout  a single  Jet.  In  Fig.  4 
we  have  made  a plot  of  tb~t£  vs.  r,  at  tb  for  all  of  the  copper  lined 
charges  studied  here  (Charge  Nos.  1-8).  Also,  in  Fig.  4,  we  have 
plotted  the  data  for  the  81.3mm  BRL  Standard  Charge  and  the  105mm  un- 
confined charge  studied  in  [1].  Three  data  points  were  selected  from 
each  jet.  With  the  exception  of  one  charge,  the  38.1mm,  90°  copper 
(Charge  No. 4),  all  of  the  data  points  for  tb~tf  fall  within  a narrow 
region  monotonically  increasing  with  increasing  r, . We  may  speculate 
that  the  90*  charge  does  not  fall  in  this  region  Because  the  low 
collapse  velocities  in  a charge  with  a large  cone  angle  results  in 
lower  pressures  and  velocities  at  the  stagnation  point  during  forma- 
tion. Therefore,  the  properties  of  this  jet  may  be  different  from 
those  of  smaller  angle  cones.  Figure  4 indicates  a correlation  be- 
tween the  two  quantities  for  the  other  charges  however.  In  fact,  if 
we  draw  a line  through  these  points,  we  obtain  a useful  breakup  time 
vs.  radius  curve,  which  can  be  used  for  design  purposes.  For  example, 
suppose  we  design  a new  copper  lined  charge  and  wish  to  obtain  its 
breakup  time  distribution.  We  can  compute  the  radius  of  the  elements 
of  the  jet  as  functions  of  time.  If  we  then  plot  the  radius  vs.  time 
curve  of  a particular  jet  element  on  the  breakup  time  vs.  radius  co- 
ordinates, we  find  that  this  curve  will  intersect  the  breakup  curve. 
This  intersection  gives  the  breakup  time  of  that  element.  We  then 
repeat  this  procedure  for  a series  of  jet  elements.  This  then  yields 
a breakup  time  distribution  for  the  new  charge.  This  breakup  time  in- 
formation may  then  be  used  together  with  other  information  on  a pene- 
tration analysis  to  evaluate  the  final  performance  of  the  new  design. 


* 


These  plots  for  the  remaining  10  charges  may  be  found  in  [9]  and  [10]. 


JET  BREAKUP  TIME 


& & 


Figure  4.  Breakup  time  vs.  jet  segment  radius  for 

various  charges.  Ten  charges  are  included. 
Three  data  points  selected  from  each  charge. 
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III.  STABILITY  OF  SHAPED  CHARGE  JETS  AS  A CAUSE  FOR  BREAKUP 


A.  Background 

The  concept  of  liquid  jet  stability  in  classical  hydrodynamics 
has  been  studied  by  many  authors  (e.g.  [11}-[14J).  The  idea  of  apply- 
ing a stability  approach  to  shaped  charge  jets  was  motivated  by  the 
remarkable  resemblance  of  the  breakup  of  shaped  charge  jets  to  that  of 
liquid  jets.  This  resemblance  is  shown  in  Fig.  5 where  the  breakup  of 
a glycerine-water  jet  [15]  is  contrasted  with  a typical  shaped  charge 
jet.  Further,  we  note  that  the  average  i/d  ratio  for  all  of  the  copper 
jet  segments  studied,  4.86  (see  Table  II),  is  very  close  to  the  value 
of  the  critical  X/2tq  ratio  of  4.5  predicted  by  classical  stability 
analyses  of  liquid  jets,  where  X is  the  wavelength  of  the  surface  dis- 
turbance, and  Tq  is  the  radius  of  the  undisturbed  jet.  These  facts 
suggest  that  the  shaped  charge  breakup  phenomenon  may  be  caused  by 
surface  instability. 

The  first  stability  analysis  of  a liquid  jet  was  published  in 
1879  by  Lord  Rayleigh  [11].  He  considered  an  inviscid  fluid  Jet  mov- 
ing at  a constant  uniform  axial  velocity  subjected  to  surface  distur- 
bances about  the  equilibrium  position.  He  obtained  the  most  unstable 
wavelength  and  the  perturbation  growth  rate  using  energy  considera- 
tions. His  values  of  wavelength  and  growth  rate  have  been  verified  by 
experiments,  and  his  work  is  still  considered  the  foundation  for  the 
study  of  jet  stability. 

Weber  [12]  extended  Rayleigh’s  work  to  include  viscous  Newtonian 
fluids.  The  wavelength  of  the  most  unstable  surface  disturbance  for  a 
viscous  fluid  does  not  deviate  too  much  from  the  ideal  fluid  case  of 
(\/2tq)  =■  4.5.  He  also  found  that  Newtonian  viscosity  tends  to  dampen 
Che  instability.  Anno  [16,17]  gave  a more  general  derivation  of  the 
analyses  by  Rayleigh  and  Weber. 

Goldin,  et  al  [15]  studied  the  breakup  of  non-Newtonian  viscous 
fluid  jets.  They  used  essentially  the  same  approach  as  Levich  [13] 
and  obtained  results  on  critical  wavelength,  growth  rate,  and  breakup 
time  for  fluids  possessing  general  viscoelastic  stress-3train  re- 
lations. In  [15],  it  is  shown  that  the  Newtonian  fluid  is  the  most 
stable  one  among  viscous  fluids.  Experimental  evidence  for  different 
types  of  fluids  seems  to  verify  their  results. 

In  our  initial  work,  we  have  determined  that  three  effects  may 
cause  jet  breakup:  surface  tension,  aerodynamic  forces,  and  strength 
of  the  jet  material.  We  have  applied  all  the  essential  results  of 
the  classical  analyses  to  the  conditions  of  a typical  shaped  charge 
jet.  We  have  also  examined  many  experimental  breakup  records  in  the 
light  of  stability  considerations.  These  results  have  ruled  out  the 
importance  of  the  surface  tension  and  aerodynamic  effects  in  the 
stability  of  a shaped  charge  jet.  This  leaves  the  effect  of  jet 
strength.  Since  strength  effects  are  not  readily  amenable  to  analyt- 
ical treatment,  we  have  turned  to  numerical  approaches  to  study  the 
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LIQUID  JET  SHApQ3  CHARGE  JET 

75  % glycerine-25  % water  copper 

REFERENCE  15  81.3MM  BRL  STANDARD 

(RADIOGRAPH  COURTESY  R.  JAMESON,  BRL) 


FIGURE  5.  The  resemblance  between  the  breakup  of 
a liquid  jet  and  a shaped  charge  jet. 
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effect*  of  jet  strength. 


In  the  next  section  of  this  report,  we  will  first  present  results 
of  s preliminary  analytical  study.  Then  more  detailed  results  on  the 
numerical  study  will  be  presented,  including  the  study  of  various  sur- 
face disturbances,  time  of  disturbance  initiation,  strain  rate,  and 
inertia  effects  in  Jets  with  strength. 

B.  Preliminary  Analytical  Stability  Study  of  Shaped  Charge  Jets 

The  results  and  formulas  developed  in  the  classical  analyses  of 
liquid  jet  stability  were  applied  to  the  case  of  a typical  shaped 
charge  copper  jet.  Analyses  for  non-stretching  and  stretching  jets 
were  used.  The  results  of  this  study  are  summarized  in  Table  IV  and 
will  be  briefly  discussed  below. 

1.  Non-stretching  Jets 

(a)  Surface  Tension.  In  our  initial  analysis,  we  have 
applied  the  classical  formulas  of  Rayleigh  (11]  and  Weber  (12]  ..o  the 
conditions  of  a typical  copper  shaped  charge  jet.  Rayleigh's  work 
governs  the  case  of  an  ideal  fluid  under  the  effect  of  surface  tension 
and  Weber's  results  are  applicable  to  a viscous  fluid  under  the  effect 
of  surface  tension.  The  growth  rate  of  disturbances  on  the  surface 

of  the  jet  as  predicted  from  these  formulas  is  quite  small  as  compared 
to  that  observed  in  radiographs  of  shaped  charge  Jets.  In  fact,  accord- 
ing to  these  formulas  the  amplitude  of  the  Initial  disturbance  only 
grows  8%  In  100  usee,  whereas  it  Is  observed  from  experiments  that  the 
shaped  charge  jet  used  for  this  analysis  breaks  up  at  approximately 
100  usee.  The  value  of  surface  tension  for  copper  used  was  1.0  N/m. 

(b)  Aerodynamic  Force.  We  have  applied  the  formulas  of 
Levlch  [13]  to  the  problem  of  air  passing  over  a shaped  charge  jet.  A 
verv  large  disturbance  growth  rate  was  predicted  by  this  analysis. 

We  feel,  however,  that  the  analysis  is  not  too  realistic  for  the  pres- 
ent problem  since  Levich  only  considers  linear  incompressible  aero- 
dynamics and  the  present  case  is  actually  in  the  hypersonic  regime. 

We  would  also  like  to  indicate  two  experimental  observations  which 
demonstrate  that  aerodynamic  effects  are  not  important  in  the  breakup 
of  shaped  charge  jets.  Vitali  [18]  has  pointed  out  that  superplastic 
jets  do  not  breakup  along  their  length  in  the  typical  manner  of  a 
copper  jet,  yet  both  jets  are  subject  to  the  same  aerodynamic  force. 
Also,  Frey  [19]  has  studied  photographs  of  copper  jets  in  a vacuum, 
and  found  the  typical  surface  disturbance  growing  within  100  usee,  sim- 
ilar to  those  found  in  jets  traveling  through  air. 

2.  Stretching  Jets 

Mikarai,  et  al  (14]  has  developed  an  analytical  approach  to  study 
the  growth  of  disturbances  on  the  surface  of  a stretching  viscous 
thread  surrounded  by  another  viscous  medium.  We  have  modified  this 
analysis  to  make  it  applicable  to  the  case  of  a shaped  charge  jet.  As 
shown  in  Table  IV,  appreciable  growth  rates  and  reasonable  A/2r  ratios 
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lalytical  Approaches  to  Shaped  Charge  Jet  Stability 
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were  predicted  by  this  enalyeia  at  time*  in  the  regime  of  typical 
breakup  times:  (Note  that  in  the  stretching  jet  case  X/2r  changes 

with  time  thus  making  the  comparison  more  complicated.)  However,  the 
approach  of  [14]  has  neglected  inertia  forces  to  make  the  governing 
equations  analytically  tractable.  Because  of  the  high  strain  rates 
present  in  the  shaped  charge  problem  Inertia  effects  will  be  Important. 
This  importance  of  Inertia  effects  will  be  verified  by  numerical  cal- 
culations in  the  next  sub-section.  Thus,  we  feel  that  the  results  of 
this  analytical  method  of  [14]  as  applied  to  the  shaped  charge  problem 
are  not  conclusive. 

C.  Numerical  Study  of  Shaped  Charge  Jet  Stability 

In  the  previous  sub-section  we  have  studied  surface  tension  and 
aerodynamic  effects  on  jet  stability  by  relatively  simple  analytical 
techniques.  Since  effects  such  as  material  strength  are,  at  present, 
not  readily  amenable  to  analytical  treatment,  we  have  used  numerical 
techniques  to  study  this  and  other  effects.  One  advantage  of  numerical 
studies  over  analytical  studies  is  that  in  the  numerical  treatment  all 
of  the  pertinent  effects  may  be  included  at  the  outset,  whereas  ana- 
lytical treatments  necessitate  the  use  of  certain  simplifying  assump- 
tions. 


1.  General  Approach 

The  numerical  study  of  jet  stability  was  undertaken  through  the 
use  of  the  two-dimensional  code  HEMP  [20,21].  The  HEMP  code  is  a 
general  purpose  code  which  solves  the  conservation  equations  of  two- 
dimensional  elastic-plastic  flow  in  plane  coordinates  or  in  axisyra- 
tnetric  coordinates.  The  solution  is  by  the  method  of  finite  differ- 
ences and  uses  the  Lagrangian  formulation.  The  code  has  the  capabil- 
ity of  handling  many  various  boundary  and  initial  conditions. 

Karpp  [22]  first  applied  HEMP  code  calculations  to  the  problem  of 
a stretching  elastic-plastic  jet  with  the  surface  slightly  disturbed. 
After  calculating  various  wavelengths  of  the  surface  disturbance,  he 
found  a broad  range  of  most  unstable  wavelengths  with  reasonable 
growth  rates.  We  have  applied  the  same  basic  method  to  study  the 
breakup  problem  in  more  detail.  The  boundary  and  initial  conditions 
cf  this  method  will  now  be  briefly  described. 

A stretching  shaped  charge  jet  is  modelled  by  a prismatic  circu- 
lar bar  fixed  at  one  end,  with  the  other  end  moving  at  a constant 
velocity.  A linear  velocity  distribution  in  the  axial  direction  is 
imposed  as  the  initial  condition,  and  the  surface  of  the  bar  is  ini- 
tially perturbed  in  the  shape  of  a cosine  function,  as  shown  in  Fig.  6. 
The  perturbed  surfaces  are  free  from  any  tractions  and  the  end  sur- 
faces are  free  of  any  shear  stresses.  Let  the  axial  velocity  be 
V (x, r, t) , the  radial  velocity  be  u(x,r,t),  the  stress  vector  on  the 
lateral  surface  be  ^(x,rg,t),  and  the  shear  stresses  on  the  end  sur- 
faces be  Trx(0,r,t),  t xg ( 0 , r , t ) , Trx(L,r,t),  Tx0(L,r,t),  then  the 
boundary  conditions  are 
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(1) 


V(Q,r,t)  - 0 

V (L, r, t)  - V0 
£<x,r,,t)  - £ 

Trx(0,r,t)  - xx0(O,r,t)  - T^CL.r.t)  - T^CL.r.t)  - 0 

and  Che  Initial  conditions  are 

V(x, r,0)  - (V0/L)x 

u(x, r ,0)  - 0 

The  calculation  may  often  be  limited  to  only  one  cycle  of  the  surface 
wave  along  the  axial  direction  because  of  the  symmetry  of  the  problem.* ** 

2.  The  Effect  of  Yield  Strength 

In  the  analytical  studies  discussed  previously  all  of  the  jets 
were  assumed  to  be  fluid  in  nature  and  the  driving  forces  for  the  in- 
stability were  restricted  to  surface  tensions,  viscosity,  and  aero- 
dynamic forces.  Now  we  will  study  the  effect  of  yield  strength  on  the 
stability  of  an  elastic-plastic  jet. 

To  study  this,  we  computed  four  jet  segments  having  different 
yield  strengths  but  being  otherwise  identical.  In  these  calculations 
we  selected  a copper  jet  segment  with  one  cycle  surface  disturbance 
having  an  initial  length  of  7.5mm  and  an  initial  mean  radius  of  1.5mm. 

The  initial  amplitude  of  the  wave  was  taken  as  10%  of  the  initial 
radius.  The  velocity  difference  between  the  two  ends  of  the  segment 
was  0.0625  mm/usec.  In  a real  jet  this  initial  time  corresponds 
typically  to  a time  50  usee  after  the  segment  was  formed  or  approxi- 
mately 100  usee  after  the  detonation  of  a BRL  81.3mm  charge.  This 
initial  configuration  and  calculation  mesh  is  shown  in  Fig.  7.  In  all 
cases  the  copper  equation  of  state  as  specified  in  HEMP  was  used. 

The  value  of  density  used  was  8.9x10^  kg/m^  (8.9  g/cm®)  and  the  elas- 
tic shear  modulus  was  taken  as  4.56x10^"  Pa  (456  kbar) . The  first  of 
the  four  segments  has  no  strength,  i.e.  purely  fluid;  the  next  three 
have  yield  strengths  of  2x10^  Pa  (.2  kbar),  2x10®  Pa  (2  kbar),  and 
2x10^  Pa  (20  kbar).  The  configuration  after  32  usee  for  each  of  these 
four  cases  is  also  shown  in  Fig.  7.  To  compare  results  we  have  com- 
puted a quantity  u>  , the  relative  growth  of  amplitude,  which  is  defined 
as  the  difference  between  the  amplitude  of  the  disturbance  after 
stretching.  A,  and  the  initial  amplitude,  Ag,  divided  by  the  initial 
amplitude,  i.e.  (A-Aq)/Aq.  This  quantity  is  also  given  for  each  case 
in  Fig.  7.  We  observe  from  these  calculations  that  the  disturbance 
grows  faster  for  jets  with  larger  values  of  yield  strength. 

* We  are  studying  a typical  element  in  the  middle  of  the  jet.  This  element  is 
being  stretched  by  the  force  of  the  forward  part  of  the  jet  on  one  side  and 
Che  rear  part  of  the  jet  on  the  other.  The  details  of  how  the  momentum  and 
energy  is  transferred  from  one  portion  of  the  jet  to  another  are  currently 
being  studied. 

**  Note  that  Karpp  and  Simon  [23]  have  recently  studied  the  strength  in 
shaped  charge  jets  using  the  experimental  results  of  rotating  charges 
and  also  by  using  two-dimensional  numerical  calculations. 
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Y * 2 x IO#Pa 
A = + 194  % 

Y = 2 x |0*Pa 
A = +247% 

t = 32/isec 

INITIAL  length:  7.5mm 

INIT  IAL  MEAN  RADIUS:  1 . 5MM 
VELOCITY  AT  MOVING  END:  0 . 0525mm/  fJ.  SEC 
A=  RELATIVE  GROWTH  OF  DISTURBANCE  = (A-Aq)/ao 
Y = YIELD  STRESS  ( EL AST  I C'PERFECT  PLASTIC) 

FICLP.i:  7.  Ci'np.i  r i si  n of  st  rt-trh  ins;  jots 
with  various  vivid  strengths. 


For  static  applications,  we  are  used  to  the  concept  that,  when 
comparing  two  materials,  the  one  with  higher  strength  will  sustain 
larger  stress,  larger  strain,  and  more  stretching.  Our  present  re- 
sults indicate,  however,  that  under  dynamic  conditions,  the  opposite 
case,  which  seems  contrary  to  common  sense,  prevails.  That  is,  under 
dynamic  conditions  of  stretching,  the  material  with  higher  yield 
stress  is  more  unstable,  will  neck,  more,  and  break  sooner. 

3 . The  Effect  of  Disturbance  Initiation  Time 

In  this  sub-section  the  effects  of  initiating  the  disturbance 
at  various  times  in  a realistic  shaped  charge  jet  are  presented.  Con- 
sider an  element  of  a jet  at  several  stages  from  the  time  when  it  is 
first  formed  until  the  time  of  its  breakup  as  shown  schematically  in 
Fig.  8.  This  element  of  jet  has  an  initial  velocity  difference  of 
0.13  mm/ysec  between  its  two  ends.  The  initial  length  and  velocity 
gradient  were  selected  such  that,  at  the  breakup  time,  the  element 
will  have  a length  equivalent  to  the  average  measured  length  of  a 
typical  jet  segment  from  a 81.3mm  BRL  standard  charge.  We  may  there- 
fore consider  only  one  wavelength  disturbance  over  this  selected  ele- 
ment. When  this  jet  element  is  first  formed  one-dimensional  calcula- 
tions indicate  its  length  whould  be  approximately  1.1mm  and  its  radius 
should  be  approximately  3mm.  The  times  chosen  to  introduce  the  dis- 
turbance, labelled  tj  through  t5  in  Fig.  8,  correspond  to  times  0 usee, 
7.92  usee,  17.62  usee,  23.69  nsec,  and  45.71  usee  after  the  element 
was  first  formed,  respectively.  The  overall  dimensions  of  this  ele- 
ment at  these  different  times  are  summarized  in  Table  V. 


Table  V 


limensions  of  Jet 

Element  Studied  at 

Various 

Times  After 

Its  Formation 

Run 

No. 

Time  after 
formation 
(ysec) 

Element 

Length 

(mm) 

Element 

radius 

(mm) 

A/2r 

Amplitude  of  Strain  Rate 

Imposed 

Disturbance  (usee  ) 

(mm) 

1 

0 

1.10 

2.98 

0.18 

0.149 

0.118 

2 

7.92 

2.13 

2.13 

0.5 

0.107 

0.061 

3 

17.62 

3.39 

1.69 

1.0 

0.085 

0.038 

4 

25.69 

4.44 

1.48 

1.5 

0.074 

0.029 

5 

45.77 

7.05 

1.17 

3.0 

0.059 

0.018 
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TRAJECTORY 
OF  THE  BACKj 
OF  THE  / 
ELEMENT  h 


LENGTH  OF  BROKEN 
SEGMENT 


-TRAJECTORY 
OF  THE  FRONT 
OF  THE  ELEMENT 


•DISTURBANCE  INITIATED 
IN  THE  ELEMENT  AT  EACH 
OF  THESE  TIMES 


-ELEMENT  FIRST  FORMED 


POSITION  IN  THE  JET 


FIGURE  8.  Schematic  showing  a position-time  plot 
of  a jet  element  which  stretches  to  the 
typical  breakup  length.  Disturbances 
are  initiated  seperately  at  each  of  the 
times  indicated  for  independent  HEMP 
calculations  of  the  element. 


The  HEMP  calculations  were  performed  using  each  of  these  times  as  a 
separate  starting  point.  The  disturbance  was  Introduced  independently 
at  each  of  these  times  and  each  was  treated  as  a separate  problem. 

The  equation  of  state  and  other  constants  for  copper  stated  in  the 
previous  problem  were  again  used.  The  yield  stress  was  2x10^  Pa  (2  kbar). 
The  amplitude  of  the  imposed  disturbance  was  taken  at  5X  of  the 
radius  in  each  instance  and  is  also  given  in  Table  V. 

The  five  HEMP  calculations  Indicate  that,  during  short  times 
after  the  initiation  of  the  disturbance,  the  amplitude  of  the  surface 
disturbance  in  Run  No. 4 grows  the  largest  amount.  This  can  be  seen  in 
Fig.  9,  where  the  relative  growth.  A,  at  a time  18  usee  after  the  dis- 
turbance initiation  for  each  case  is  plotted  versus  the  initial  strain 
rate.  This  gives  an  indication  of  how  soon  each  disturbance  begins  to 
grow.  As  time  proceeds,  however,  and  we  continue  the  calculations, 
the  earlier  configurations  (Run  Nos.  1-3)  eventually  reach  a con- 
figuration similar  to  Run  No. 4,  and  then  the  disturbance  begins  to 
grow  quickly.  This  can  be  seen  in  Fig.  10  where  the  relative  growth 
of  Run  Nos.  2,3,  and  4 are  plotted  versus  time.  Examining  Run  No. 2, 
we  see  that  the  disturbance  is  stable  until  a time  of  approximately 
30-35  usee  after  formation,  then  the  growth  rate  drastically  in- 
creases. Run  Nos.  3 and  4,  in  which  the  disturbance  is  initiated  at 
times  later  than  Run  No. 2,  also  begin  to  grow  more  rapidly  in  the  re- 
gion 30-35  usee.  Thus  we  see  that,  no  matter  how  early  we  initiate 
the  disturbance,  when  a critical  time  is  reached,  the  disturbance  will 
begin  to  grow  and  eventually  the  same  final  breakup  configuration  will 
be  obtained. 


4.  The  Effect  of  Disturbance  Wavelength 

The  effect  of  disturbance  wavelength  was  studied  by  consider- 
ing a jet  of  realistic  radius  and  strain  rate  at  a reasonably  early 
time  after  formation  and  introducing  surface  disturbances  of  differ- 
ent wavelengths.  A jet  element  having  an  initial  nadius  of  2.13mm 
and  an  initial  strain  rate  of  0.061  usec-i  was  used.  The  equation  of 
state  and  material  constants  for  copper  stated  in  the  previous  problem 
were  again  used.  Four  separate  surface  disturbances  having  initial 
wavelengths  of  1.065mm,  2.13mm,  4.26mm,  and  8.52mm  were  introduced  and 
HEMP  calculations  made  for  each  case.  Note  that  we  have  chosen  this 
problem  such  that  the  case  X » 2.13mm  will  eventually  grow  to  the 
observed  segment  length  after  the  experimentally  determined  breakup 
time  for  a typical  shaped  charge  jet.  Let  us  denote  this  "correct" 
initial  wavelength  as  Xq  » 2.13mm.  We  can  then  denote  the  other 
cases  as  i * Xq/2,  * = 2.' g and  X * 4*0'  The  relative  growth,  as 
computed  using  the  HEMP  code  for  each  of  these  cases,  is  plotted  vs. 
time  in  Fig.  11.  From  this  plot  we  see  that  the  disturbances  having 
initial  wavelengths  of  X * *0/2  and  X » 4Xq  grow  very  slowly.  The 
wavelength  X = Xq,  which  will  grow  into  the  proDer  segment  size,  grows 
very  quickly  but  not  quite  as  fast  as  the  case  X = 2Xg. 
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RUN  NO. 4 


INITIAL  STRAIN  RATE  (fiSEC  ' ) 


FIGURE  9.  Relative  growth  vs.  initial  strain 
rate  for  five  successive  element 
configurations  at  18usec  after 
disturbance  Initiation  in  each 
conf igurat ion. 


26 


5.  Irregular  Surface  Disturbance 

To  examine  these  wavelength  effects  further,  we  have  numerically 
studied  J.ets  having  an  Irregular  surface  disturbance.  In  previous 
cases,  we  have  considered  Jets  having  a surface  of  a single  sinusoidal 
curve.  Next,  we  have  combined  the  four  wavelengths  of  the  previous 
study  to  obtain  an  irregular  surface  disturbance  which  is  described 
by  the  following  function: 
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where  k^  ■ Xg/2,  k2  * Xg,  k3  “ 2Xg,  ^4  ■ 4Xg.  We  have  taken  numer- 
ical values  to  be  rg  ■ 2.13mm,  Xg  - 2.13mm,  Ag  « 0.026mm,  and  the 
initial  strain  rate  equal  to  0.061  usec“l.  We  again  used  the  material 
constants  for  copper  described  previously.  The  initial  configuration 
which  has  a length  of  2Xq  is  shown  in  Fig.  12,  together  with  the  con- 
figuration after  52  usee.  We  observe  that  the  jet  appears  to  neck  in 
two  places.  This  indicates  that  Xg  is  the  dominant  or  fastest  grow- 
ing wavelength,  which  is  in  agreement  with  the  experimentally  ob- 
tained segment  length.  To  demonstrate  this  more  quantitatively,  we 
have  fit  a Fourier  series  to  the  outer  surface  cf  the  configuration 
after  52  usee  shown  in  Fig.  12  and  have  examined  the  coefficients  of 
the  terms  of  each  wavelength.  Denoting  the  length  after  52  usee  as 
2Xj,  we  observe  that  the  largest  of  all  the  Fourier  coefficients  is 
that  of  the  term  containing  the  length  ■ • ■ A comparison  of  the 
coefficients  of  the  four  wavelengths  of  interest  are  also  shown  in 
Fig.  12.  Thus , this  result  indicates  that  the  component  with  an 
initial  wavelength  of  X » v0*  which  stretches  into  * ■ X^  after 
52  usee  is  the  most  critical  one. 


In  the  case  above,  toe  initial  disturbance  consists  of  four 
waves  of  different  wave  1 er.  gt  h , which  a re  .1 11  "in  phase".  We  performed 
yet  another  irregular  surface  disturoance  calculation.  In  this  calcu- 
lation, a "random"  disturbance  was  used  which  was  comprised  of  five 
waves  of  different  wavelength  and  phase  angle--,  so  that  the  waves 
would  not  be  "in  phase"  at  the  ends.  Thus,  this  new  disturbance 
represents  more  closely  a random  one.  The  function  used  for  the  sur- 
face was 
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irregular  surf  act1  disturbance. 


FIGURE  13.  Results  of  HEMP  calculations  of  a stretching  jet 
with  "random"  surface  disturbances. 


6.  The  Effect  of  Inertia 

To  study  the  importance  of  inertia  effects  on  the  growth  of 
surface  disturbances  in  stretching  shaped  charge  jets,  we  have  made 
HEMP  code  calculations  of  three  jets.  Each  jet  was  identical  in  all 
respects  except  for  the  density  of  the  material  used  in  the  calcu- 
lations. The  configuration  used  was  the  standard  single  wave  surface 
disturbance  depicted  in  Fig.  6 with  rg  - 2.13mm  and  wavelength  equal 
to  2.13am.  The  Initial  amplitude  used  was  Aq  - 0.107mm  and  the 
Initial  strain  rate  equal  to  .061  usee-*.  The  equation  of  state  and 
material  constants  described  previously  for  copper  were  used  except 

were  used  pq  » 1x10*  kg/m*, 
kg/m-*.  The  relative  growth  vs. 
time  for  each  case  is  plotted  in  Fig.  14.  We  observe  that  increased 
density,  i.e.  increased  inertia  force,  has  a retarding  effect  on  dis- 
turbance growth. 
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now  three  different  values  of  density 
pq  ■ 8.9x10*  kg/m*,  and  pq  ■ 16.5x10* 


Flf'.UKF.  14.  Relative  growth  vs.  time  for  the  stretching  jots  of 
different  densities  is  predicted  by  HEMP. 


IV 


CONCLUSIONS 


The  major  results  and  conclusions  of  the  present  research  may 
be  summarized  as  follows: 

A.  Measurement  of  Jet  Radiographs 

1.  Data  from  jet  radiographs,  including  jet  velocity, 
breakup  time,  and  segment  size,  for  jets  from  eleven 
charges  were  compiled  and  plotted  for  easy  reference. 

2.  These  measured  jet  velocity  and  jet  radius  data 
agree  closely  with  those  calculated  from  a one-dimensional 
shaped  charge  model. 

3.  Analysis  of  this  data  together  with  using  one- 
dimensional  calculations  resulted  in  a semi-empirical 
breakup  time  V9.  jet  radius  curve  for  copper  jets.  This 
curve  can  be  used  to  estimate  an  approximate  breakup 
time  for  new  charge  designs. 

B.  Stability  of  Shaped  Charge  Jets 

1.  Simple  stability  analyses  and  experimental  results 
indicate  that  surface  tension  and  aerodynamic  forces 
are  not  important  in  the  breakup  of  shaped  charge  jets. 

2.  Numerical  HEMP  code  calculations  of  stretching 
elastic-plastic  jets  subjected  to  surface  disturbances 
were  conducted.  The  effects  of  material  yield  strength, 
time  of  disturbance  initiation,  wavelength  of  the  distur- 
bance, irregular  surface  disturbances  and  jet  density 
were  examined.  The  following  results  were  found: 

a.  Jets  with  higher  yield  strengths  break  sooner,  all 
else  being  equal. 

b.  Jets  with  lower  densities  will  break  sooner,  all 
else  being  equal. 

c.  For  the  shaped  charge  jet  calculated,  there  is  a 
critical  time  for  the  growth  of  the  disturbance 
amplitude.  Disturbances  introduced  early  do  not  grow 
appreciably  before  this  time,  but  grow  rapidly  after 
this  critical  time.  Disturbances  introduced  after  this 
time  also  grow  rapidly. 

d.  A critical  wavelength  (or  a range  of  wavelengths) 
exists;  disturbances  having  this  wavelength  grow  faster 
than  all  others.  The  length  of  the  broken  jet  segment 
caused  by  this  critical  wave  is  in  the  range  of  measured 
jet  segment  lengths. 

e.  When  irregular,  or  random,  surface  disturbances  are 
introduced,  the  growth  is  again  dominated  by  the  distur- 
bance component  with  the  critical  wavelength.  The  jet 
surface  grows  into  a shape  similar  to  that  obtained  if 
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only  Che  wave  with  Che  critical  length  were  Introduced. 

3.  Materials  with  high  density  and  low  yield  strength 
are  promising  liner  materials  which  will  likely  retard 
the  breakup  of  the  resulting  jet. 
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Appendix  A 


This  eppendix  contains  Che  specifications  of  the  eleven  charges 
studied  in  this  report. 
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FIGURE  A2.  Geometry  of  the  50.8  mm,  42°  Copper  lined 
charges  with  different  wall  thicknesses. 
(Charge  Nos.  5-7) 

(Courtesy  of  R.  .Jameson) 
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Appendix  B 


This  appendix  contains  a graphical  display  of  jet  velocity,  Jet 
breakup  time,  and  Jet  radius  data  computed  from  the  jet  radiographs 
of  the  eleven  charges  studied  in  this  report. 
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FIGURE  B2.  Theoretical  and  experinental  jet  velocity  vs. 

jet  particle  position  for  t h»-  series  of  38.  liam, 
wall,  a 1 umi nun  iined  charges  with 
various  cone  angles,  at  various  times  after  ttie 
arrival  of  the  detonation  wave  at  the  cone  apex. 
(Charge  Nos.  9-11) 
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position  tor  the  50.8mm,  *2  , 1.524mm  wallfeopper  line 
charge  at  130.5  usee  after  detonation  wave  arrival  at 
apex.  (Charge  No.  6) 


FIGURE  US.  lhoorotic.il  and  experimental  jet  velocity  vs.  jet 

particle  position  for  the  50.8mm,  42  , 2.54mm  wall, 
copper  lined  charge  at  130.7  psec  after  detonation 
wave  arrival  at  cone  apex,  (charge  No.  7) 


HCURE  Bf>.  Thooret  ical  and  exper  iment  a 1 jet  velocity  vs.  Jet  particle  position  for 
tin-  81.3mm,  43°,  tapered  wall, Copper  lined  charge  at  130.5  nsec  after 
detonation  wave  arrival  at  cone  apex.  (charge  No.  8) 
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FICURE  B8.  Jet  breakup  time  vs.  jet  particle 
position  at  t»93.4  usee  for  the 
38.1mm,  40  , l.lfiS  mm  wall  copper 
lined  charge. 
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FIGURE  B9.  Jet  breakup  time  vs.  jet  particle  position 
at  t = 1 2 2 . 7 psec  for  the  38.1nm,  80°, 

1.168  mm  wall  Copper  lined  charge. 

(Charge  No.  3) 


/ 


300  400  500  600  700  800 

JET  PARTICLE  POSITION  £ (mm) 


Jet  breakup  time  vs.  jet  particle  position 
at  t*133.3  psec  for  the  38.1cm,  90  , 
1.168  mm  wall  copper  lined  charge. 

(Charge  No.  4) 
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FIGURE  BIO. 
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FIGURE  Bll.  Jet  breakup  time  vs.  jet  particle  | 

position  at  t=U3.1  usee  for  | 

the  38.1mm,  40°,  1.626mm  wall  aluminum  ^ 

lined  charge. 

(Charge  No. 9) 
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FIGURE  B12.  Jet  breakup  time  vs.  jet  particle  position 
at  t=97.7  psec  for  the  38.1mir,  60  , 
1.626  mm  wall  aluminum  lined  charge. 
(Charge  No.  10) 
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FIGURE  B13.  Jet  breakup  time  vs.  jet  particle  position 
at  t= 1 32 . 8 usee  for  the  38.1mm,  90  , 
1.626  mm  wall  aluminum  lined  charge, 
(charge  No.  11) 
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FIGURE  B16.  Jet  breakup  time  vs.  jet  position  at 
t«130.7  usee  for  the  50.8  mm,  42°,  2.54  ram  wall, 
copper  lined  charge  (charge  no.  7). 
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FIGURE  B18.  Theoretical  and  experimental  jet  radius  vs. 
jet  particle  position  at  t«107.7  usee  for  the  38.1  mm,  20°, 
1.168  naa  wall,  copper  lined  charge  (charge  no.  1). 
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FIGURE  B19.  Theoretical  and  experimental  jet  radius  vs.° 
jet  particle  position  at  t=93.4  Msec  for  the  38.1  mm,  40°, 
1.168  mm  wall,  copper  lined  charge  (charge  no.  2). 
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FIGURE  B21.  Theoretical  and  experimental  jet  radius  vs.  jet 
particle  position  at  t»133.3  -sec  for  the  38.1  mm,  90°, 

1.168  mm  wall  copper  lined  charge  (charge  no.  4). 
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FIGURE  B22.  Theoretical  and  experimental  jet  radius  vs.  jet 
particle  position  at  t“113.1  usee  for  the  38.1  mm,  40°, 

1.626  mm  wall  aluminum  lined  charge  (charge  no.  9). 
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FIGURE  B23.  Theoretical  and  experimental  jet  radius^vs.  jet 
particle  position  at  t-97.7  usee  for  the  38.1  mm,  60°, 

1.626  mm  wall  aluminum  lined  charge  (charge  no.  10). 
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FIGURE  B24.  Theoretical  and  experimental  jet  radius  vs.  jet 
particle  position  at  t"132.8  ysec  for  the  38.1  mm,  90°, 

1.626  mm  wall  aluminum  lined  charge  (charge  no.  11). 
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FIGURE  325.  Theoretical  and  experimental  jet  radius  vs.  jet 
particle  position  at  t=130.6  psec  for  the  50.8  mm,  42°, 

0.762  mm  wall  copper  lined  charge  (charge  no.  5). 


FIGURE  B26.  Theoretical  and  experimental  jet  radius  vs.  jet 
particle  position  at  t=130.5  psec  for  the  50.8  mm,  42°, 

1.524  mm  wall  copper  lined  charge  (charge  no.  6). 


JET  RADIUS  (m> 


PARTICLE  POSITION  IN  JET  «(m) 


FIGURE  B27.  Theoretical  and  experimental  jet  radius  vs.  jet 
particle  position  at  t=130.7  usee  for  the  50.8  mm, 

42°,  2.54  mm  wall  copper  lined  charge  (charge  no.  7). 


PARTICLE  POSH  1(11  III  THE  JET  CW 


FIGURE  B28.  Theoretical  and  experimental  jet  radius  vs. 
jet  particle  position  at  t=130.5  Msec  for  the  81.3  mm 
42°,  tapered  wall  copper  lined  charge  (charge  no.  8). 

(Note:  The  lead  particle  at  1003  mm  in  the  jet  has  a 
radius  of  4.6  mm  and  is  not  plotted  here  for  convenience.) 


Appendix  C 


This  appendix  contains  tabulated  data  for  the  jet  velocity,  break- 
up time,  length,  diameter,  radius,  aspect  ratio  and  AVj  of  the  jet 
segments  as  computed  from  radiograph  measurements  for  the  eleven 
charges  studied  in  this  report. 

The  breakup  times  in  these  tables  are  referred  to  t=0  when  the 
charge  is  first  initiated.  Throughout  the  rest  of  the  report,  including 
the  plots  of  appendix  B,  all  times  are  referred  to  t=0  when  the 
detonation  wave  first  reaches  the  apex  of  the  liner  cone.  The  time 
it  takes  for  the  wave  to  travel  from  the  point  of  initiation  to  the 
apex  of  the  cone  for  charges  5,6  and  7 is  6.4ysec;  for  all  other 
charges  this  time  is  7.0ysec. 
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